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ABSTRACT 

We review the observational evidence for extra mixing in stars on the red giant 
branch (RGB) and discuss why thermohahne mixing is a strong candidate mechanism. 
We recah the simple phenomenological description of thermohaline mixing, and aspects 
of mixing in stars in general. We use observations of M3 to constrain the form of 
the thermohaline diffusion coefficient and any associated free parameters. This is done 
by matching [C/Fe] and [N/Fe] along the RGB of M3. After taking into account a 
presumed initial primordial bimodality of [C/Fe] in the CN-weak and CN-strong stars 
our thermohaline mixing models can explain the full spread of [C/Fe]. Thermohaline 
mixing can produce a significant change in [N/Fe] as a function of absolute magnitude 
on the RGB for initially CN-weak stars, but not for initially CN-strong stars, which 
have so much nitrogen to begin with that any extra mixing does not significantly affect 
the surface nitrogen composition. 

Subject headings: (Galaxy:) globular clusters: individual (M3), stars: abundances, 
stars: evolution, stars: Population II 



1. Introduction 

Standard stellar evolution theory predicts that only one mixing event will change the surface 
composition of a low -mass star as it ascends the red giant branch. That event is the so-called first 



dredge-up (FDU, see llbenlll967l ) associated with the inwards migration of the base of the convective 
envelope into regions where hydrogen burning via the CNO-bicycle has occurred. Relatively modest 
changes in surface C and N abundance are predicted, and once the convective envelope recedes 
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outwards these changes are brought to a halt. However, obse rvations of low-mass red giants (M < 
2.5M0, see for example ICharbonnel &: Do Nascimentdll998l ) for essentially all compositions show 
trends among light element abundances that cannot be accounted for by the FDU. Some form of 
non-convective mixing seems to occur whereby greater amounts of the products of partial hydrogen 
burning are cycled into the convective envelope over a much longer timescale, and during more 
advanced phases of RGB evolution, than can be explained by the FDU. The observational results 
summarised below imply that current canonical models do not include essential physics of this 
so-called "extra mixing" . Regardless of what the physical mechanism is, extra mixing is required 
to conform to the following observational criteria: 



1. It commences after the hydrogen burning shell has erased a composition discontinuity in the 



event, and may continue to at least the tip of the RGB ( 


Gilrov & BrowJl991 


Charbonnel et al. 


1998; 


Gratton et al. 


2000: 


Smith & Martell 


2003: 


Shetrone 2003: 


Weiss & Charbonnel 


2004 





Martell et al.ll2008l ). The onset of the extra mixing is thus thought to coincide with a lo- 
cal maximum (the so-called "bump") observed in the RGB luminosity function of globular 
clusters. 



It must occur over a range of masses and metallicities (jSmiljanic et al.ll2009l . and references 
therein), being activ e in giants of all metallicities fr om solar to at least j Fe/H] ~ —2.5 
(jGratton et al.ll2000l ) and masses less than ~ 2.5Mq (jLambert fc Rieslll97?l ). although not 
necessarily with equal efficiency throughout these mass and metallicity ranges. 



It must deplete ^Li (jCharbonnel et al 



1998 



Smilianic et al.ll2009l : iLind et al.ll2009l ). 



It must decrease the ^^C/^^C ratio ( Charbonneill994 . 1996), since values lower than predicted 



by the FDU are found amon g Population I field g iants ( 



19811: 



200c 



Charbonnel et al.lll998 '). open cluster giants ( Gilrov 1989 : Gilrov fc Brownlll991 : Smilianic et al 



20071 ) and halo field giants (ISneden et al.lll98()l : iGratton et al.ll2000l ). 



bmkin et al 



1976; 



Lambert &: Ries 



MiKolaiti. .1 j2nidrdobular clust er fct,.t,ronmR.dn-RI„,n.o T.v.rnJ 



5. It must decrease the total carbon abundance since systematic decreases with advancing lu- 



halo field giants ( 


Suntzefi 


1981 




1989 


: Carbon et al. 


1982 




Trefzeer et al. 


1983: 


Laneer 


et al. 


1986: 


Gratton et al.l 


2000 




Bellman et al. 


2001 




Smith & Martell 


2OO3I: 


Martell et al. 


20081: 


Shetrone et al. 


20ld) 


. In Population II g 


jiants the behaviour of the carbon abundance can 



serve as an even more potent probe of the extent of extra mixing than the ^^C/^^C isotope 
ratio, because the latter can attain near-equilibrium values for only moderat e amounts of 
mixin g that would otherwise cause only small (~0.1 dex) changes in [C/H] (jSneden et al 



19861 ). 



6. 



As a consequence of the previous point it must increase the nitrogen abundance. The results 
of CN cycling are observed on the upper half of the red giant branch. Halo field stars on the 
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upper RGB were found by iGratton et al.l ([20001) to show an excess of nitrogen compared to 
those on the lower RGB. 



It is expected tha t the mechanisni(s) w i ll also destroy ^He inside the star (Dearborn et al 
198^ iHata et al.l Il995l : iDearborn et al.l Il996l : ISackmann &: Boothrovdl Il999l : ICharbonnel &: ZahnI 



2007al lbh. As we cannot observe ^He in stellar atmospheres directly this is not an observational 
constraint, but it is a significant requirement from the study of chemical yields and galactic evolu- 
tion. The importance of ^He is discussed in Section 3. 

In the last four decades many candidate extra-mi x ing mechanisms have been suggested. Thes e 



include: rotatio nal mixing (ISweigart &: Mengell Il979l: 



magnetic fields (jPalmerini et a. 



20091 : 

ity waves ( Denissenkov &: Toutll20od ). Individually, none of these have been proven to be satisfac- 
tory. Due to its promising ability to account for the above requirements, and the necessity of its 
occurrence in low mass giants just after the FDU, in this study we focus on ^He-d riven "thermoha- 
line mixing" (Eggleton et al-lbood . EDL06 hereafter, Charbonnel fc Zahn 2007a . CZ07a hereafter, 
Eggleton et al.ll2008l . EDL08 The name comes from a phenomenon seen in oceans, and 

is taken from the two major determinants of the density of sea water - its heat content ("thermo") 
and its salinity (the salt or "haline" content). It is common to find warm salty water overlying 
cool fresh water. Although the higher salinity of the warm water makes it denser, the higher heat 
content acts to stabilise the stratification. The subsequent evolution of the system is determined by 
the competition between two diffusion processes and their associated timescales - the time for the 
(stabilizing) heat to diffuse away compared to the time for the (destabilizing) salt to do the same. 
Hence the process is often called "doubly-diffusive" , and it has been studied in the oceanographic 
context for many years (see recent reviews, theoretical modelli ng, observations and labor atory ex- 
periments in Progress of Oceanography Volume 56, 2003; e.g. iRuddick fc: GargettI l2003l ) . Within 
oceans it is now well known that the rapid diffusion of heat from the warm salty layer produces 
an over-dense layer that begins to sink into the cooler fresh water below. The temperature stays 
roughly the same as the surrounds, and "salt fingers" form which extend downward delivering the 
saltier water to deeper regions. Reciprocal fresh-water fingers move upward and replace the salty 
water with fresher water. Laboratory experiments have also played a role in helping to characterise 
the instability . Wor k by Stommel and Faller published in ISternI (Il960l ) as well as more recently 
Krishnamurtil (|2003l . and references there in) have helped elucidate the instability. 



Chaname et al 



Nordhaus et al. 



20081 : 



2005: 



Busso et al 



PalacioseLall l2006|) 



2003), and internal grav- 



A similar process can occur in stars. He re it is not sa l t but the mean molecular weight that 
is the "destabilising agent", in the words of iDenissenkovi (120101 ). Usually, the molecular weight 
increases as we move toward the centre of the star, as a result of nuclear burning and fusion 



^One must concede that multiple processes an d indeed interactions between them can affect the transpo rt. Models 
have been made that include multiple processes ( Cantiello fc Langer 20ld : Charbonnel fc LagardellioiB ) by simply 
adding the diffusion coefficients for each proces s. This does not allow for the interaction between the processes as 
discussed bv lDenissenkov &: Pinsonneault 
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reactions. If it were to decrease, then the plasma would be buoyantly unstable, just as is the case 
in convection. However, just as in the oceanic case, we must include the rapid thermal diffusion 
which can act to stabilise the motion. In this paper we consider the situation where some local 
event causes a decrease in the molecular weight in an otherwise stable region within a star. We 
investigate the effects of the resultant "therm ohaline mixi i ig" or doubly diffusive process that is 
initiated by a molecular weight (/i) inversion (jUlrichI Il972l : iKippenhahn. Ruschenplatt &: Thomas 
Il98' 



Re cently thermohaline mixing has featured prominently i n the literature. EDL06, CZ07a, 
EDL08, ICantiello k. Langeii (|20ld ) and lCharbonnel fc Lagardd (j20ld ) have discussed in detail the 
important consequences of its inclusion during the RGB. The dichotomy in RGB carbon abun- 
dance between metal poor stars and carbon-enhan ced metal-poor star s has b een explai n ed by 



Stancliffe et al 



(I2OO9I') using therm ohaline mixing. ICantiello fc Langerl (j201Cl ). IStancliffd (I2OIOI ) 



and ICharbonnel &: Lagardd (j2010|) have shown that its operation beyond the giant branch may 
affect the subsequent asymptotic giant branch (AGB) evolution. This mechanism may be a crucial 
part of stellar physics that has been missing from the models. We thus believe it pertinent to 
investigate the effects of the mechanism in some detail. 

In this paper we will take the approach of trying to model the mixing in spherically symmetric 
stellar models and investigate its ef fect on observable surface ab undances. In this respect our 
approach follows that of CZ07a and ICharbonnel &: Lagardd (120101 ). We will examine the change 
of the abundances of carbon and nitrogen on the red giant branch of globular clusters, focusing on 
the case of M3. 

Currently a range of approaches is taken to include thermohaline mixing in evolution codes, 
especially if the ^^C/^^C ratio is the constraint used to determine the extent of extra mixing. The 
12(-j^13q j-g^^^jQ -^yg^g Qj^g Qf ^]^g ^Tsi indicato rs that extra mixing must operate on the RGB. It is 



1975 



Tomkin et al 



1976 



Charbonnel 



class ically used to probe the results of FDU (jPearborn et al 
1994 ). It naturally follows that the ^^C/^^C ratio could be used to trace the extent of extra mixing 



and constrain any mechanism. The change in ^^C/^^C ratio following FDU will depend on the 
efficiency of mixing and allow us to explore the mixing velocity via a diffusion approximation. 
EDL08 estimated the mixing speed with their formula for the diffusion coefficient and found that 
a range of three orders of niagnitude in their free parameter can lead to the observed lev els of 
and ^He depletion, luirichi (|l972l ) and Ixippenhahn. Ruschenplatt fc Thomas! (|l98d ) both 
use essentially the same formula (UKRT formula hereafter) for the diffusion coefficient but their 
choice of the free parameter varies by two orders of magnitude. In an attempt to constrain the 



^This has been referred to as mixing" by EDL06 and EDL08, to distinguish it from a separate occurrence 
of "thermohaUne mixing" in stars. In that case we may have mass transfer in a binary system, where material of 
a higher mean molecular weight is accreted on an envelope of lower molecular weight. This situation is unstable 
and some thermohaline circulation will take place to redistribute the composition of the star to result in a stable 
stratification (exchanging energy with the thermal content of the m aterial in doing so). Thi s case is not relevant t o 
the discussions in this paper, and more information may be found in lChen fc Hani (|200J) and lStancliffe et al.l (|2007|). 
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parameter space we address the following questions: 



1. The ^^C/^^C ratio is generally used as a tracer to probe the extent of mixing. This ratio 
saturates near the CN equilibrium value rather quickly in low metallicity stars, and hence is 
of limited utility. Is there a better way to constrain the mixing? 

2. Which formalism should be used? Here we will limit our investigation to the EDL08 and 
UKRT prescriptions for the diffusion coefficient. 

3. Once the preferred formalism is identified, what diffusion co-efficient (or mixing velocity) is 
needed to match observations? What values do we use for any free parameters? 



As has been t he practice fo r many years, we turn to globul ar clusters to test ou r understanding 
of stellar theory. ISmithI ^QM), ISmith k MartelJ (|2003l ) and iMartell et al.l (|2008l ) have compiled 
observations of carbon and nitrogen along the giant branch of M3. This has provided us with a 
valuable alternative to the ^^C/^'^C probe. By matching our models to the carbon depletion (as 
a function of absolute magnitude) observed in this cluster we can attempt to constrain both the 
form of the thermohaline diffusion coefficient and the values of any parameters contained therein. 
As carbon and nitrogen are intrinsically linked in the CN burning cycle we include observations of 
nitrogen as an additional tracer. Furthermore, we identify when extra mixing begins in the models 
and compare this to the observed luminosity function bump (LF bump) in the cluster. 



2. Thermohaline Mixing 

The usual condition for convective instability is simply that a blob moved from its equilibrium 
position will be buoyantly unstable and continue to move away from its initial position. In a 
region of homogeneous chemical composition this results in the usual Schwarzschild criterion for 
instability: 

Vrad > Vad, (1) 

where Vad = (91nr/51nP)ad is the adiabatic gradient and Vrad = (51nT/(?lnP)i.ad is the same 
gradient assuming all the energy is carried by radiation (one usually uses the diffusion approximation 
for this expression). This equation simply says that the steepest stable gradient is the adiabatic 
one, and if a steeper gradient is required to carry the energy by radiation, then radiation will fail 
and buoyancy will develop. 

This condition however ignores the possibility of a variation in the chemical composition. 
Including this leads to the Ledoux criterion for instability: 



Vrad > Vad + (|) V^, 



(2) 
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where (p = (91np/91n fi)p^T, ^ = —{dlnp/dlnT) and = {din fi/dlnP). The essential feature 
here is the appearance of the term involving the gradient of the molecular weight (the multiplicative 
terms (j) and 6 come from the thermodynamics and are of little importance for the present discussion; 
both are unity for a perfect gas without radiation pressure). 

For thermohaline mixing we need the Ledoux criterion to be broken, but with the added 
condition that the molecular weight gradient decreases with depth, i.e. 



< 0. 



(3) 



Thermohaline mixing has a long history. IStothers fc SimonI (jl969l ) were the first to consider 
how thermohaline mixing may impact upon the structure and subsequent evolution of a star. They 
theorised that accreted material from a helium-rich companion would be thermohaline unstable 
and may be responsible for the pulsations of f3 Cepheids. 



Soon after [Abraham &: IbenI (j 19701 ) and lUlrichI (jl97ll ) recognised an unusual property of the 
reaction 

(4) 



^He(^He,2p) He, 



which is that despite being a fusion reaction, it actually reduces the mean molecular weight because 
it produces three particles from two. (These particles also have more kinetic energy than the 
initial two particles, because the reaction is exothermic.) They thought the reaction may have an 
important role during pre-main-sequence contraction, and could lead to the situation where < 
and a thermohaline instability may develop. This was later found to have little effect due to the 
short time scale of the pre-main-sequence. 



UlrichI (I1972I ) was the first to derive an expression for the turbulent diffusivity in a perfect gas. 



He considered thermohaline mixing during the cor e flash where off-centre ig nition leads to carbon- 
rich material sitting on material that is helium-rich. iKippenhahn et al.l (|l98d ) extended this to allow 
for a non-perfect gas which included radiation pressure and degeneracy. Thermohaline mixing is 
unlikely to h ave an appreciable effe ct during the core flash because of the very short timescale 
of the flash. iDearborn et al.l (|2006l ) showed that a small amount of overshooting inwards could 
remove the molecular weight inversion on a dynamical timescale, thus removing the possibility of 
any thermohaline mixing. 



3. Thermohaline Mixing on the First Giant Branch 

A further application, and the one of interest to us here, was proposed by EDL06. During 
main-sequence evolution a low mass star produces a substantial amount of ^He. Close to the centre 
all of the ^He has been destroyed in completing the pp chains. But at lower temperatures the 
^He remains (e.g. see Fig 1 of EDL08). When the FDU occurs it mixes the envelope abundances 
over a large region, producing a homogeneous envelope with a ^He content that is far above the 
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equilibrium value. Following FDU, the convective envelope recedes, leaving behind a region that is 
homogeneous in composition. 

When the hydrogen burning shell advances, the fragile ^He begins to burn and the resulting 
decrease in the molecular weight leads to an inversion developing. Note that ^He burning alone is 
not sufficient to drive the mixing. Whenever there is hydrogen burning via the pp chains the ^He 
+ ^He reaction (Equation will decrease fi. However the combined effect of the other reactions 
and the low abundance of ^He makes the effect of the ^He + ^He reaction negligible. If one can 
homogenise the region beforehand, then the effect of ^He + ^He dominates. This way thermohaline 
mixing is initiated at essentially the same luminosity as the bump in the luminosity function, since 
both are caused by the H-shell reaching the abundance discontinuity left behind from FDU. 

From equation (4) it can be shown that the change in fi must be (EDL08) 

6fi = fi^6{X{^Be))/6. (5) 

and they found the magnitude of such inversions to be of the order A///// ~ 10~^. Although this 
seems small, convection is in fact driven by a similarly small superadiabaticity. 

We now try to understand what happens to the material in this region of the star. When ^He 
burns, a parcel forms that is slightly hotter and has lower molecular weight than its surroundings. 
This means it has a higher pressure than it requires for its position. Hence it quickly expands (and 
begins to cool) in order to establish pressure equilibrium. The expansion reduces the density and 
therefore the element becomes buoyant. The parcel rises until it finds an equilibrium point where 
the external pressure and density are equal to that inside the bubble. This is expected to be a 
small displacement which occurs on a dynamical timescale. 

As the molecular weight inside the bubble is lower than its surroundings the equilibrium point 
must correspond to a place where the external temperature is higher than that of the bubble. The 
temperature inside the bubble will be lower than its surroundings (we may assume a perfect gas 
equation of state): 



where subscript i denotes the inside of the bubble and subscript o denotes the surroundings. As 
heat begins to diffuse into the parcel, we expect layers will start to strip off in the form of long 
fingers. It is this secondary mixing that governs the overall mixing timescale. The mixing cycles in 
fresh ^He from the envelope reservoir, to replace the ^He-depleted material that is rising toward the 
star's convective envelope. This upward-flowing material will also have experienced other burning, 
such as CN cycling, while in the hotter region. It will also experience further burning in the future 
as the material is cycled through this region on subsequent occasions. 

The diffusion of heat and composition is analogous to the situation in the oceans where haline 
fingers diffuse salt and heat, although in the stellar context we are dealing with a compressible 
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flo w. For a tho r ough investigation of the difi'erences between the oceanic and stehar case we refer 



Denissenkovl (|20ld ). The fact that in both instances there are two diff^usive processes at work 



to 

has seen the term 'thermohahne' adopted in the astrophysics hterature (e.g. CZ07a). 



As stated above, EDL08 originaUy referred to thermohahne mixing on the RGEij as "(5/i mixing" 
to distinguish and emphasise the mechanism that drives the mixing, namely the slight difference 
in the molecular weight brought about by '^He burning. Their simulations seemed to hint at the 
possibility of a dynamical phase in addition to the linear phase which CZ07a modelled as UKRT 
thermohahne mixing. 

Many groups have realised the importance of understand ing the mixing and are invest i gating 
the instability that arises from the ^He burning. The anlysis by lDenissenkov Sz Pinsonneaultl (120081 ) 
of the mixing sug gests a dynamical p hase should not arise. Further complicating the picture are 
the 2D models by iDenissenkovl (|2010l ) that predict a different behaviour in the fluid compared to 
the ID parameter fitting of CZ07a. 

We consider 5fj, mixing to be a strong extra-mixing candidate along the RGB as it can 
meet the criteria outlined in the introduction (CZ07a, EDL08). To summarise, this ^He burn- 
ing drives mixing that begins after the hydrogen shell encounters the composition discontinuity 



cleosynthesis and Big Bang Nucleosynthesis remain consistent (Dearborn et al. 


198d 


Hata et al. 


1995; Dearborn et al. 1996; Sackmann & Boothrovd 1999: Charbonnel & Zahnl2 


007a|t 


). Measure- 


ments of ^He in HII regions and planetary nebula (Rood et al. 1984 


Hoean 1995; 


Charbonnel 1995: 


Charbonnel & Do Nascimento 1998; Tosi 1998: Palla et al. 200ol: ] 


ilomano et al 


2003 


) match the 


predicted yields from Big Bang Nucleosynthesis. Canonical models predict that low- 


mass, main- 



sequen ce stars are net producers of ^He which is returned to the ISM through mass loss. iHata et al 
(jl995l ) have shown that about 90% of the ^He produced on the main sequence must be destroyed 



to reconcile the two fields. We saw above how these stars produced ^He, and the newly discovered 
mechanism now also destroys it, removing the inconsistency with Big Bang Nucleosynthesis. This 
will still allow for the exi s tence of the minority of planetary nebula that show large ^He abundances 



(|Balser et al.lll997l . Il999l . 



2006 



Charbonnel fc Zahnll2007bl V 



The thermohahne mixing mechanis m described here will operate until the ^He is destroyed 



which may be beyond the tip of the RGB (|stancliffJ201ol ; ICantiello fc Langerlboiol :! ICharbonnel fc Lagardd 



2010l ). The molecular weight inversion caused by the ^He burning creates an instability that cy- 
cles CN processed material into the envelope, having the des i red e ffect on the surface abundances 
(CZ07a). Finally CZ07a, EDL08 and ICharbonnel Lagardd (j2010l ). demonstrated that the degree 
of mixing depends on mass and metallicity, most clearly seen in the different values of the carbon 
isotopic ratio in Population I and Population II stars. 



It s applic ation to the AGB is s till co ntentious, see iKarakas et all (|201ol ). Istancliffel l|20ld) . ICantiello &: Langen 
( 20101 ) and ICharbonnel fc Lagardd l|201ol ). 
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4. The Formula for the Diffusion Coefficient 



It is common in stellar interior calculations to use a diffusion equation to simulate mixing. It 
is important to remember that many mixing processes, such as convective mixing, are advective in 
nature, not diffusive. In the former, a property is distributed as a result of bulk flows within the 
fluid. The latter follows from Pick's first law, which postulates that a flux exists between regions 
of high and low concentration of the quantity of interest. The equation governing this is 

j = -DVn (7) 

where n is the density of the quantity and j is its flux. The parameter D is the "diffusion coefficient" 
of the "diffusivity" . For particles with a mean-free-path / and typical speed v then 

D = ^vl. (8) 

As no theory for time dependent mixing exists, it is common to use a diffusion equation to simulate 
mixing in stellar interiors. Historically, this is also how thermohaline mixing has been calculated. A 
diffusion equation is solved for each species in the star, and thus determines the radial composition 
variation. 

EDL08 used the following formula, based on estimates of the mixing velocity and the convective 
formalism in their evolution code: 



D 



(/i - /imin) if {k > k^in) 

(9) 



if {k < /Cmin), 

where k is the mesh point number, counted outwards from the centre of the model, A;min is the 
value of k for which /i takes its minimum value of Umim f is the radial coordinate, F\av is a constant 
which is selected to obtain the desired mixing efficiency and tnuciear is an estimate of the nuclear 
evolution timescale (see EDL08). This is an extension of the convective formula used in the EDL08 
calculations, with the addition of a — /Xmin term to reflect the driving by the inversion. Thus 
this formula is not a local one as the value of D at a point in the star depends on the value of /imin 
at some other location in the star. Nevertheless, it is a phenomenological form that recognises the 
role of the /i inversion in driving the mixing. 

This formulation ensured the correct region was mixed but note that the mixing speed is 
formally zero at the position where ^ has its minimum even though the mixing should presumably 
be the most efficient at this point. EDL08 give upper and lower estimates for the mixing velocity 
and find that they can change their free parameter -Finy, and consequently alter the speed by 
three orders of magnitude whilst still producing the observed levels of ^^C/^^C and leading to ^He 
depletion. 
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CZ07a adopt the UKRT formulation for thermohaline mixing, which results from a linear 
analysis of the problem. They cast it in the following way, 

where K is the thermal diffusivity and Ct is a dimensionless free parameter. In fact Ct is related 
to the aspect ratio, a, of the fingers in the following manner: 



Ct = -vr a 



2„2 



(11) 



The appropriate value to use for this parameter remains uncertain. To understand thermoha- 
line mixing fully clearly requires a hydrodynamic theory that will determine the diffusion coefficient 
and any associated parameters. In lieu of such a theory, we can compare with observations and see 
what form of D and values of constants are needed to match the real world. This is the approach 
we adopt in this paper. 

Alternatively, one could try to determine Ct or a by comparison to laboratory experiments 
and the oceanic case, which has been well studied. Stomme l and Faller carried out experiments 
of fluids in laboratory conditions. This work, published in ISternI (119601 ) . saw the devel opment 
of lon g salt fingers with lengths that were larger than their diameters (a ~5). This led lUlrich 
(|l972l ) to determine that Ct — 658 which was the basis of the choice by CZ07a to use Ct = 
1000. The huge differences between incompressible salty water and the th e self-gravitating pl asma 
we are considering make direct comparisons very difficult, however (see iDenissenkovl l20ld ). In 
contrast to the oceanic case, Kippenhahn envisaged the classical picture where mixing is due to 
blobs. If Kippenhahn 's expression is cast into the same fori n as Equation 1101 t hen Ct ~ 12. This 
value seems to agre e with the 2D-hydrodynamical models of IDenissenkovl ()2010l ) and ID models of 
Cantiello &: Langen (j2010l ). but is inconsistent with the value preferred by CZ07a and the present 
work (see below). 



Models and Results 



Messier 3 was chosen as a test case for thermohaline mixing in Population II stars because 
i n many respects it can be considered a typi cal globular cluster. A metallicity of [Fe/H] = — 1.4 



( Sneden et al. 



2004 



Cohen &: Melendeal2005l ) means that it falls near the mode in the metallicity 



distribution of halo globular clusters. Among clusters of similar metallicity, the ho rizontal branch 
morphology of M3 ha s stars on both the red and blue sides of the RR Lyrae gap (jSandagelll970l : 
Buonanno et al.lll994l ). and the colour distribution is intermediate be tween extreni e s such as those 



of M 13 and NGC 6752 (whose HBs have more extended blue tails; iNewelll Il970l : iLee &: Cannon 



1980l) and NGC 7006 (whose HE stars are more concentrated to the red side of the RR Lyrae 



gap; 



Sandage &: Wildeylll967l ). Red giants with enhanced A3883 CN bands in their spectra were 
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discovered in M3 by Suntzeff (1981), and further studied by iNorris fc SmithI (jl984l ). The CN 



distribution contains both CN-strong and CN-weak giants (i.e., high and low nit rogen abundance 



giants) as is typical of other clusters of similar metallicity (jNorris &: SmithI Il984l ). The behaviour 



of the CN bands at a given lumi n osity on the RGB appears to anticorrelate with carbon and 



oxygen abundance (jSuntzefll Il98ll : ISmith et al.l 119961 ) . When giants of different luminosity are 



compared a significant trend becomes apparent in the carbon abundance. Among red giants with 



( Suntzefi 


1981, 


1989: 


Smith 


2002) 



field giants ([Smith &: Martelll l2003l ) . It is this behaviour of carbon in particular that is of great 
use in constraining models of thermohali ne mixing. The carbon isotope ratio ^^C/^^C has a low 
value of ~ 4-6 among the bright giants (jPilachowski et al.l l2003l : iPavlenko et al.l |2003| ) , although 
the number of stars for which this has been measured is small. In addition, M3 displays O, Na, 
and Al abundance inhomogeneities of t he type that are commonplace within globular clusters of 
similar metallicity (jCarretta et al.ll2005l). The cluster has a population of stars that are enhanced 



in Na and Al but depleted in oxygen (iKraft et al.lll992l : ICavallo fc Nagarj|2000l : ISneden et al. 



2004 



Cohen fc Melendezlbood : I.Tohnson et ahlbood ^. but without the extremes ir i O abundance depl etion 



and Na enhancement that are found in the oft-compared cluster M13 (jSneden et al.ll2004l ). A 
substantial component of the 0-Na-Al variations in MS are arguably the products of early cluster 
self-enrichment, and thermohaline mixing on the first as cent giant branch is not expec ted to have 
played a role in their production. Work by CZ07a and ICharbonnel &: Lagardd (|20ld ) also show 
that these elements are not affected by thermohaline mixing. 



We use MONSTAR (The Monash version of the Mt.Stromlo evolution code; see lCampbell fc Lattanzio 



20081 ) to produce stellar models for M3. The stars are evolved from the zero-age main-sequence un- 
til the core flash. They are of mass M=0.8 Mq, metallicity Z=0.0005 and have a solar-scaled CNO 
abundance, as a first approximatio n. Populati on II subdwarfs in th e globular cluster metallicity 
range tend to have [C/Fe] ~ (e.g. lLairdll985l ). In M13 the work of iBrilev et al.l (|200J) shows an 
upper envelope to the data that approaches [C/Fe] to —0.2 near the main sequence. 

Our mixing length parameter a is set to 1.75 and we run without any overshoot. Mixing is 
calculated using a diffusion equation, and we investigate different formulae for the diffusion co- 
efficient used for thermohaline mixing. These parameters approximate the s tars in MS quite well : 
MS has an age estimated by various authors to be between 11 .S to 14.2 Gyr (jChabover et al.l Il992l : 
Jimenez et al]ll996l : IvMidenBer j[200ol : ISalaris Weisslbood : IaIvcs et alTl2004l l and [Fe/H]=-1.5. 
Figure [1] shows the evolution of our model star in the HR Diagram. We mark the end of core H 
exhaustion (CHE in the figure, at an age of t=12.5 Gyr), the beginning and end of FDU (bFDU and 
eFDU respectively) as well as the position of the LF bump. Figure [2] shows a Kippenhahn diagram 
for the evolution of our star. Here we see the convective envelope moving inwards (coinciding 
with the beginning of FDU), whilst the deepest point of penetration marks the end of FDU. The 
location of thermohaline mixing is the lightly shaded region between the envelope and the advancing 
hydrogen burning shell. In Figure [Sj using the UKRT formulation for thermohaline mixing and 
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Ct=1000, we plot the variation of the surface [C/Fe] as well as the ^^C/^^C ratio versus My. Here 
it can clearly be seen that the ^^C/^^C ratio drops below 10 very quickly, finally reaching a s low 
as 6.5. This is in go od agreement with the observations of 4-6 by iPilachowski et al.l (j2003l ) and 
Pavlenko et al.l (|2003l ). The gradual depletion of carbon throughout the entire RGB on the other 
hand provides a more sensitive constraint on the mixing. 

Figure m compares the results with the EDL08 and UKRT formulae for the diffusion coefficient. 
In this figure open circles denote CN-weak stars, filled circles denote CN-strong stars and crosses 
represent CN-intermediate stars all from M^. In this figure as well as the proceeding figures the 
errors in magnitude are smaller than the symbols used. The ran dom errors in the abundance 
measurements can be up to ±0.3 dex according to ISuntzefa (|l98ll ). although this does not take 
into account potential systematic errors due to limitations in the stellar atmospheres code used in 
the abundance derivations. We have run four models to compare to the observations, two with 
thermohaline mixing that use the UKRT formula for the diffusion coefficient and two models that 
use the EDL08 formula. The dashed line and the solid l ine are variations of the U KRT prescription, 
the broken line uses Ci=12 (a ~ 1) as suggested by iKippenhahn et al.l (Il980l ) while the model 
represented by the solid line uses Ct=1000 (a ~ 6), a geometry more akin to CZ07a. The dotted 
and dot-dashed lines (which nearly sit on top of each other) are EDL08 models where i<inv=100 
corresponds to the dotted line and Finv=10 corresponds to the dot-dashed line. The two choices of 
Finv lie within the three orders of magnitude that lead to the required level of ^^C/^^C depletion 
as discussed in EDL08. 

The two EDL08 models destroy the available ^He very quickly and are unable to cycle much 
CN processed material to the envelope. As previously mentioned, in the EDL08 formulation of 
the diffusion coefficient the mixing is dependent on the difference in /i, i.e. (// — /imin)- Material 
near the envelope is cycled down faster and the mixing becomes less efficient closer to the location 
where /i has its minimum. Material is not replenished at the same rate it would be with a local 
condition on the /i gradient as seen with the UKRT prescription. This in turn affects the depth at 
which the minimum value of /i occurs. We can see in Figure [5] that the location at which fj, has its 
minimum occurs further out in the EDL08 scheme. The steep temperature gradient ensures the 
fragile '^He is easily destroyed while the carbon-rich material is not exposed to temperatures where 
it can burn significantly. CN cycling is therefore reduced whilst the driving fuel, ^He, is easily 
processed irrespective of one's choice of -Finv The bottom panel in Figure [5] clearly demonstrates 
the dependence of the location of the minimum value of /i on the mixing scheme. The four models 
in Figure H] are plotted with the same symbols in Figure [5j We have plotted the velocities, carbon. 



*We note tha t ther e appears to be a lack of CN-strong stars on the lower RGB. This is an artifice of the 
oriRinal Suntzefj (198l|) study in which the lower luminosity stars that were observed happened to be CN-weak. 



Norris fc SmlthNigS^lh showed that CN-strong giants do exist in M3 at luminosities corresponding to the faint limit 
of the Suntzeflf i 1981 \ survey. However, their later study concentrated on the A3883 CN band strength and did not 
measure the [C/Fe] or [N/Fe] abundances. Consequently, absence of CN-strong stars at the faintest limits of the 



abundance plots in our paper is a consequence of observational effects. 
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nitrogen and molecular weight profiles at a hydrogen-exhausted core mass of Mc=0. 32970 Mq 
corresponding to My = —0.34 . The bottom panel demonstrates the depth to which the models 
mix. This highlights the fact that the two UKRT models have their // minimum occurring closer 
to the hydrogen shell and mix deeper than the EDL08 models at the same core mass. As expected 
the two EDL08 abundance profiles show little CN processing. As shown by EDL08 this is sufficient 
to produce the observed ^^C/^^C values but as Figure [5] shows it is not enough to reproduce the 
[C/Fe] or [N/Fe] variations. Note that to calculate a velocity from Equation [8] we require a value for 
I. We have used / = 1.75x Hp where Hp is the pressure scale height. This may not be an accurate 
estimate but it is expeted to give us an indicative velocity. 

The surface abundances are determined by both the depth to which the material is mixed 
as well as the speed of mixing. The depth of mixing in a "steady state" is determined by the 
location of the /i minimum which is in turn dependent on the mixing speed, so that the two are 
not independent. Consider the case of very fast mixing, such as convection. The abundance would 
be homogeneous from the deepest point of mixing to the surface. Clearly the surface abundance is 
controlled by the burning conditions at the bottom of the mixed region - a case we call "burning 
limited." Now consider the case of very slow mixing. The abundance profile is marginally altered 
from the radiative case. Although the material is mixed to the surface it proceeds so slowly that 
the advance of the burning region is almost independent of the (slow) mixing - just as it would be 
in the limiting case of no mixing (i.e. a radiative zone). Hence in the slow mixing case the surface 
abundances are "transport limited." In Figure H] the bottom panel suggests that the UKRT Ct=12 
model is transport limited and this, not the depth to which material is mixed, is responsible for the 
lack of processing. The bottom panel in Figure [5] shows that the EDL08 models do not mix as deep 
as the UKRT models. Yet the EDL08 models destroy essentially all of the ^He (see Figured]) while 
this is not true for the UKRT formulation. ^He and consequently carbon in the Cj=12 model are 
not destroyed because not enough material is physically transported through the radiative region 
and processed before the end of the RGB. In Figure [6] we see that there is is almost a difference 
of two orders of magnitude between the mixing speeds in the two UKRT models. In addition at 
this stage of the evolution the Ct=12 model does not connect the region of the /U-inversion (where 
the burning occurs) with the envelope to cycle in fresh fuel (it does so further up the RGB): this 
further hinders any processing. 

It appears that the requirement for significant carbon depletion is two fold: the position of the 
minimum in fi must occur at a reasonably high temperature and we need a sufficiently fast mixing 
velocity. As a result of our comparison with M3 we prefer the UKRT formalism over that used by 
EDL08. 



6. Constraining the Diffusion Coefficient 



As previously mentioned iKippenhahn et al.l (Il980l ). lUlrichI (j 19721 ) and CZ07a differ in their 
choice of a, the aspect ratio of the fingers by a factor of five and this translates to a difference of 
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two orders of magnitude in Cj. We can also determine a velocity for the mixing through a diffusion 
equation. 

Having selected our preferred formalism we can now use observations to constrain the diffusion 
coefficient. In Figure Eh- we illustrate the effect of standard evolution on the carbon and nitrogen 
abundances. In all the Figure [7] panels, open circles denote CN-weak stars , filled circles denote 
CN-strong stars and crosses represent CN-intermediate stars as before. The uncertainty in the data 
is as specified in Section 5. Figures [7)3, ^ and [TJl include thermohaline mixing with the UKRT 
prescription. In each p anel we provide four mo dels, these represent different values of Ct- We 



include: Ct= 12 as per iKippenhahn et al.l (|l98Cl ) (dashed lines, a ~ 1), Ct =1000 as per CZ07a 



(solid lines, a ~ 6) and two intermediate values, (7^=120 (dot-dashed lines, a ~ 2) and Ci=600 
(dotted lines, comparable to Ulrich 1972 Ci=658, a ~ 5). 

The canonical evolution in Figure [7^ shows very little carbon depletion following FDU. This is 
contrary to the observations for the upper part of the RGB. We include this panel to highlight to 
the reader the need for extra mixing. In Figure [Tja we have assumed a scaled solar CNO abundance 
and varied the free parameter as de scribed above. It is unlikely that M3 possesses a scaled solar 



abundance. iRussell Sz Dopital (|1992| ) have already shown that Large Magellanic Cloud abundances 
are not scaled Solar and we would not expect the environment during globular cluster formation 
to resemble the solar neighbourhood. 

In Figure [7J: we have altered the initial CNO abundances whilst keeping the metallicity and 
total CNO content constant. As the carbon abundance does not change a great deal before the 
onset of extra mixing we have assumed the most carbon-rich stars are representative of the pre- 
FDU abundance, similarly for the most nitrogen-poor stars. To match the constraints, the following 
initial CNO abundances were assumed X(C)=5.45x 10-^ X(N)=1.5xlO-^ and X(0)=2.86x 10"^. 
We altered the carbon and nitrogen so that we arrived at our desired post-FDU value and adjusted 
the oxygen abundance in order to keep the total CNO constanlH. We have decreased the carbon 
by a factor of 1.47 and the nitrogen by a factor of 1.63 (in relation to solar). These values could 
easily be fine-tuned to better fit the data but it is the proof of concept we are concerned with. 
The two slow mixing cases are unable to account for the depletion in carbon. Both Cf=600 and 
Ci=1000 are a good fit to the carbon and t he nitrog en. It does appear that the mixing should 



begin at a slightly lower luminosity, however ISweigartI (|l978l ) demonstrated that lowering the He 



abundance will cause the extra mixing to begin earlier because the shell reaches the discontinuity 
sooner. We find that our models provide a fit to the upper envelope of the run of [C/Fe] vs My in 
Figure [7t but can only fit the nitrogen for the CN-weak stars. This argues that much of the spread 
in [N/Fe] in the lower panels of Figure 4 is dominated by a variation that was present among the 
cluster stars before they commenced red giant branch evolution. The fact that none of the models 
in the the lower panels of Figure 4 can reproduce the observed spread in [N/Fe] is suggestive of the 
CN inhomogeneity in M3 having a predominantly primordial origin. This inference is consistent 



^Note that the efTect of FDU is very small so the pre-FDU and post-FDU values are nearly the same. 
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with the discovery in a number of clusters, including M3 (jNorris &: Smithlll984l ). that N abundance 
variations (or equivalently CN band strength variations) are found both at and b e low the magnitude 
of the bump in the RGB luminosity function (e.g., see lSuntzeff &: Smithlll99ll : iBrilev et al.ll2004l : 
Carretta et apboosl : for the clusters NGC 6752 and M13 which have similar [Fe/H] to M3). This 
argues strongly for the spread being present in the stars at their birth. 

In Figure [7Jl we have assumed the CN-strong stars are a separate population with their own 
CNO abundance. Moreover we have assumed there is a CN-[C/Fe] anticorrelation in which the CN- 
strong stars are initially 0.4 dex depleted in carbon relative to the CN-weak stars (jSmithll2002l . and 
references therein) . The CN-strong stars show the results of hot hydrogen burning in the gas from 
which they formed, including CN and ON cycling. We have therefore changed the initial helium 
abundance from Y=0.2495 to Y=0.28 as well as the CNO valuej^. Here we set X(C)=1.9x 10~^ 
X(N)=5.5x 10~^ and X(0)=2.6x 10~^. We find in this case the slower mixing is a good fit for the 
carbon abundance whilst the faster mixing is able to account for the more extreme observations of 
carbon depletion. These models are unable to match the more extreme nitrogen enhancements in 
the CN strong stars. The [N/Fe] is initially so high even on the main sequence that any enhancement 
due to extra mixing only slightly raises the [N/Fe]. We conclude that much of the nitrogen spread 
in the CN strong stars is due to primordial variations as a result of hot hydrogen burning (in an 
early generation of stars that instigated cluster enrichment) rather than thermohaline mixing in 
the present-day giants. 

Our findings are consistent with those of CZ07a in that Ci=1000 is our preferred value of 
the free parameter in the UKRT thermohaline mixing prescription. In Figure [8] we investigate 
the effect of increasing Ct beyond 1000. The dotted line corresponds to the EDL08 prescription 
with Finv=100, the solid line corresponds to Ct=1000, the dashed line to Ct=3000 (a ~10) and 
the dotted line Ct=12000 (a ~20). Once again we use the same symbols for the M3 data. The 
initial abundances are solar scaled CNO, however it is the shape of the profiles we are interested in 
here. We include the EDL08 model to compare the similarities in the ^He depletion as a function 
of magnitude. The fact that the UKRT ^He profiles begin to approach that of the faster mixing 
EDL08 model once again suggests that depth to which the UKRT models mix is responsible for the 
carbon depletion. Initially the faster mixing cases are able to deplete the carbon more efficiently 
than our preferred value of Cj=1000. Although all three values of Ct begin to deplete at different 
rates they eventually lead to similar levels of CN processing. The similar levels of carbon depletion 
suggest that a steady state has been reached and for Ct > 1000 the process is burning limited. 



7. Mass Loss 



Mass loss could have a significant effect on the level of carbon depletion in the models. If there 
is less envelope, the amount of carbon that needs to be processed to reduce the [C/Fe] is also less. 
In Figure [9] we investigate the role of mass loss in the UKRT formula for mixing. We provide three 
models for the Ci=120 case, standard Reimers' mass loss (dotted line), no mass loss (dashed line) 
and a factor of two increase on the standard mass-loss rate (solid line). The effect of mass loss on 
the surface composition is greatest towards the tip of the giant branch where the mass-loss rate is 
highest, but even there the alteration in [C/Fe] is slight compared to a zero mass- loss model. These 
findings are consistent with EDL08, who also find the effect of mass loss negligible after the onset 
of thermohaline mixing. 



8. The Onset of Mixing and the Location of the LF Bump 



Our final observational test concerns the onset of mixing and the observed location of the 
abundance changes. The hydrogen-burning shell will reach the composition discontinuity indepen- 
dent of any thermohaline mixing, which only begins after this event. Once the mixing starts, the 
rate is dependent on the parameter Ct and so the surface abundance begins changing at a visual 
magnitude that depends slightly on Ct- The models for the four different values of Ct in Figure [7| 
commence thermohaline mixing within AMy = 0.02 mag of each other. In F igurellOlwe have taken 



Figures Eb and Eli and marked the position of M3's LF bump acc o rding to (jFusi Pecci et al 



199C 



grey dashed line) and the position according to ISmith &: Martelll (|2003l . dot-dashed grey line). In 
all four panels we also include the magnitude at which the hydroge i i shel l reaches the discontinuity 
in the CN weak model^(the s olid vertical line) . Fusi Pecci et al. ( 199C ) determine the LF bump 
of M3 to occur at My=— 0.06, ISmith Martelll (|2003l ) at My=0.45 whereas our models reach the 
composition discontinuity at My=— 0.26. In the UKRT formula for the diffusion coefficient the 
change in the surface abundances begins almost immediately aft e r the bump. These results differ 
somewhat from the non-rotating case of ICharbonnel Lagardd ([20101) who find there is a delay 
of ~ 0.5 magnitudes before the onset of any changes to the surface abundance, albeit for different 
masses and metallicities to those considered here. The differences between the two quoted values 
for the LF bump are due to the respective choices of the distance modulus. We find our models 
are closer to the older results but we refrain from drawing any conclusions. A study that includes 
lower luminosity giants is required in order for us to make further comparisons to the observations. 



In Figure [TOb we find that increasing the amount of ^He in the model has delayed the onset of 



®We note that this is consistent with the growing literature on multiple stellar populations in GC and their 
presumed *He content jpiottolbopgl . and references there in). 



^Technically, we plot t h e loca l luminosity maximum that occurs during this phase; this is Lb,n 
of ICharbonnel fc Lagarde l|2O10l l. 



in the notation 
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the bump. These results are consistent with [Sweigarg (|l978l ) who in addition to this showed that 
the size of the discontinuity is also reduced with increasing ^He abundance. This is owing to the 
fact that an increase in ^He translates to less hydrogen and less of a discontinuity. The fact that 
the hydrogen shell encounters the discontinuity later for models with higher initial ^He is due to 
the effect the composition has on the penetration of the convective envelope. A higher hydrogen 
abundance allows for deeper penetration of the envelope and hence the shell reaches this depth 
sooner. 



9. Conclusion 



M3 is a well studied system that demonstrates the abundance patterns we commonly associate 
with globular clusters. Along with many other clusters it displays significant [C/Fe] depletion along 
the RGB, the implication being that some form of internal, non-canonical mixing must be occur- 
ring. Our models with thermohaline mixing show tha t the carbon and nitrogen observations can 
be explained if we adopt the hybrid theory outlined in ISmithI (120021 ) , where stars in the cluster are 
undergoing extra mixing as they ascend the RGB and the presence of primordial abundance inho- 
mogeneities due to ON cycling are needed to explain the initial carbon and nitrogen abundances. 
We have used observations of M3 to investigate our theoretical understanding of thermohaline 
mixing. Our findings are summarised below: 



The variation of [C/Fe] with magnitude provides a much more stringent test of any pro- 
posed extra-mixing mechanism than simply matching the final ^^C/^'^C ratio. When data 
of sufficient quality is available this can constrain the details of any proposed extra-mixing 
formulation. In the present case the UKRT formulation of thermohaline mixing is a far better 
fit than the phenomenological prescription given by EDL08, although both fit the constraint 
provided by the carbon isotope ratio. 

The UKRT prescription of thermohaline mixing with C;:=1000 seems to best fit the data for 
M3. This is consistent with the results of ICharbonnel &: Lagardd (|2010l ) for higher metallici- 
ties and CZ07a for a range of metallicities. 

We infer that there is a spread of ~0.3 to 0.4 dex in [C/Fe] in the stars in M3 from their 
birth. Without this initial difference in [C/Fe] between the two populations, thermohaline 
mixing cannot reproduce the change in [C/Fe] seen on the giant branch. That there are two 
populations is absolutely required, because once Ct is sufficiently large, an increase of this 
coefficient doesn't produce a bigger A [C/Fe]. Primordial C and N inhomogeneities have been 
directly observed as abundanc e differences among main sequence star s in globular clusters 
such as M13 and NGC 6752 (|Brilev et al.l bood : ICarretta et al] bood ) . which have similar 
metallicity to M3. 
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4. Thermohalin^ mixing can produce the observed values of the carbon isotopes seen in M3. 



5. To reproduce the entire spread of [C/Fe] values seen in the giants of M3 it is essential that 
thermohaline mixing operate in both the CN-strong and CN-weak populations identified in 
(3). In this case we can explain the full sprea d in [C/Fe] seen near the tip of the giant branch 
in M3. A similar exercise was carried out by iDenissenkov &: VandenBergi (120031 ) for the case 
of M92, modelled with a simple parameterized extra-mixing formulation. The data for M92 
are from many different sources and make it difficult to estimate precisely where the extra 
mixing begins. For this reason we do not discuss M92 further in this paper. 

6. Thermohaline mixing can produce a significant change in [N/Fe] as a function of My on the 
RGB for initially CN-weak stars but not for initially CN-strong stars, which have so much N 
to begin with that any extra mixing does not significantly affect the surface composition. 

7. The level of depletion of carbon is dependent on the depth to which the material is mixed 
and how fast it is mixed. 

8. Mass loss has little effect on the surface abundances. 



9. Both the predicted an d observed composition changes take place at a luminosity that is higher 
than the LF bump in iFusi Pecci et al.l (|l990l ) . The observed abundances be gin to decrease at 
a luminosity lower than the LF bump preferred by lSmith &: Martelll ()2003l ). Uncertainties in 
the distance modulus make it difficult to draw further conclusions. 



We have seen that the linear theory of Ulrich ( 1972 ) and Kippenhahn et al. (jl98ol ) provides 
a fit to the carbon and nitrogen abundances in the giants in M3 (assuming there are two different 
populations initially). CZ07a have shown that the same theory (and indeed the same pararn eter 
Ct=1000) seems to fit field stars of a range of metallicities (see also lCharbonnel &: Lagardell2010l ). It 
is important to note that thermohaline mixing is more than a theory with an adjustable parameter. 
For example, its beginning is determined clearly by the fusion of ^He which produces a molecular 
weight inversion. Also, the hydrodynamics provides the physical formulation for the diffusion co- 
efficient used, and hence its variation throughout the star, to within a constant which depends on 
the geometry of the fingers expected in the mixing process. Nevertheless, until we have a complete 
theory that also determines this geometric factor, or at the least some numerical simulations, there 
is a gap between understanding the fundamental physics and the sort of wo rk presented here, to 
match the observations. To this end we note the work of iDenissenkovl (j2010l ) whi ch addresses this 



issue . There has been rec ently considerable progress in modelling the oceanic case ([Stellmach et al 
2010l : iTraxler et al.ll2010l ). We are ourselves working on 3D hydrodynamic models in the stellar 



context, which will be the subject of another paper. 



^ Note, when we refer to "t hermohaline mixing" we are referring to the hnear theory as proposed by Uhich 1 19721 ) 
and lKippenhahn et aL I l|l980l ). 
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Although our models can explain M3 very well the question remains whether this work can 
be extended to all clusters. The UKRT prescription for thermohaline mixing appears to model 
the internal mixing of a young metal rich cluster . Old me t al-po or clusters such as M92 will have 
undergone a different mixing history. Furthermore lSweigarti (jl978l ) suggests changing the metallicity 
an d helium conte n t will drastically alter the location and size of the LF bump. Preliminary work 
by lAngelou et al.l (120101 ) will be the subject of subsequent studies. This work further supports the 
ability of thermohaline mixing to explain extra mixing on the RGB. 
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Fig. 1. — HR Diagram for our star that models M3. We mark: core hydrogen exhaustion (CHE) 
and the beginning of FDU (bFDU) which we define as the point where the penetration of the 
envelope begins to affect the surface abundances. Also marked is the end of FDU (eFDU) and the 
bump, after which extra mixing is expected to begin. 
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Fig. 2. — A simplified Kippenhan plot for our model star. We plot the interior of the star in 
mass co-ordinates against model number which is a non-linear proxy for time. The grey shaded 
region shows the penetration of the convective envelope, its deepest point coinciding with the end 
of FDU. We also mark the top of the hydrogen-burning shell. Unshaded regions denote radiative 
zones, whilst the shaded region between the shell and the envelope represents the location where 
thermohaline mixing occurs in the star. Although not obvious on this scale there is still a small 
radiative buffer between the top of the shell and the thermohaline mixed regions. Thermohaline 
mixing has the effect of extending the convection zone and allows material to be mixed closer to 
the shell where it can be processed. 
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Fig. 3. — Top panel: ^^C/^^C as a function of magnitude. Bottom panel: [C/Fe] as a function 
of magnitude. In this figure we compare the evolution of [C/Fe] to that of ^^C/^^C. The UKRT 
thermohaline mixing formula with C(=1000 has been used in this model. It can be seen that the 
isotopic ratio reaches equilibrium soon after the onset of extra mixing, where as [C/Fe] depletes 
along the entire RGB. 
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Fig. 4. — Abundances as a function of visual magnitude. In this figure we compare the UKRT 
prescription for the diffusion coefficient to that of EDL. Open circles denote CN-weak stars in M3, 
filled circles denote CN-st r ong s tars and crosses represent stars of intermediate CN strength, with 
all data taken from lSmithI (|2002l ). The two UKRT models correspond to the dashed line for Ct=12 
and the solid line for Ct=1000. The two EDL models are denoted by the dotted line for Finv=100 
and the dot dashed line for Finv=10, these sit nearly on top of each other on this scale. The models 
are all of mass M=O.8M0 and Z=5 x 10"^ with solar scaled CNO abundances. The errors m 
magnitude are smaller than the symbols used. The maximum error in [C/Fe] and [N/Fe] is ±0.3 
dex. 
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Fig. 5. — Various profiles for the same models as Figm'e [H The same symbols are used. The 
profiles are all taken at the same stage in evolution, that is when the hydrogen-exhausted core 
mass Mc=0. 32970 Mq and at and My = —0.34. Top Panel: the mixing velocity in the region 
above the shell. Second Panel: the carbon abundance in the region above the shell. Third Panel: 
the nitrogen abundance in the region above the shell. Bottom Panel: highlights the respective 
locations where the minimum value of fi occurs. 
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Fig. 6. — Velocity, Xc and ^ profiles for the two UKRT models. This is an expanded plot of the 
first panel in Figure [5] so is once again taken at a hydrogen-exhausted core mass of Mc=0. 32970 
and My = —0.34 . The Ct=12 model (dashed line) mixes two orders of magnitude more slowly 
than the Q=1000 model (solid line). It also fails to transport material to and from the envelope 
at this stage. It does connect with the envelope before the end of the RGB. 
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Fig. 7. — Carbon and nitro gen abundanc es as a function of absolute magnitude in M3. Here we 
compare our models to the ISmithI (|2002l ) compilation of M3 observations. Open circles denote 
CN-weak stars, filled circles CN-strong stars, whilst the crosses represent stars of intermediate 
CN strength. The errors in magnitude are smaller than the symbols used. The maximum error 
in [C/Fe] and [N/Fe] is ±0.3 dex. Where extra mixing is included in our models dashed lines 
correspond to Ct=12 (a ~ 1), dot-dashed lines to Ci=120 (a ~ 2), dotted lines to Cf=600 (a ~ 5) 
and solid lines to 6^=1000 (a ~ 2). Figure [T]: includes models with the following initial CNO 
abundances X(C)=5.45x 10-^ X(N)=1.5x 10"^ and X(0)=2.86x 10"^. In Figure EH the initial 
CNO abundances are X(C)=1.9x 10"^ X(N)=5.5x 10"^ and X(0)=2.6x 10""^. 
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Fig. 8. — Abundances as function of magnitude. We consider values of Q beyond 1000. M3 data 
is represented with the same symbols as in the previous figures. The dotted line corresponds to the 
EDL prescription with Finv=100, the solid line corresponds to Ct=1000, the dashed line to Ct=3000 
(a '-^lO) and the dash-dotted Hne to Ct=12000 (a ^-^20). There is little difference in the level of 
carbon depletion once Ct > 1000, because the depletion is limited by the burning rate rather than 
the transport rate. 
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Fig. 9. — Carbon abundance relative to solar as a function of magnitude. Here we plot UKRT 
Cf=120 models with varying levels of mass loss. Symbols are described in previous figures. The 
dashed line corresponds to a model with no mass loss, the dotted line to standard Reimers' mass-loss 
rate and the solid curve twice the mass-loss rate. There is little effect on the carbon depletion. 
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Fig. 10. — Carbon and nitrogen abundances as a function of absolute magnitude in M3. These are 
the same panels as Figures [Tb and[7li with the location where the local luminosity maximum occurs 
when the hydrogen shell m eets the composi t ion di scontinuity (solid vertical line), the location of 
the LF bump for M3 from iFusi Pecci et al.l (Il990l ) (dashed vertical line) and location of the LF 
bump for M3 from lSmith fc Martelll ^200^ ) (dot-dashed vertical line). 



